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, We present the first attempt of systematically investigating the effects of shell 

OO \ correction energy for a dynamical process, which includes fusion, fusion-fission and 

quasi-fission processes. In the superheavy mass region, for the fusion process, shell 
correction energy plays a very important role and enhances the fusion probability 
^ \ when the colliding partner has a strong shell structure. By analyzing the trajectory 

in three-dimensional coordinate space with the Langevin equation, we reveal the 
. mechanism of the enhancement of the fusion probability caused by 'cold fusion 

0^ \ valleys'. The temperature dependence of shell correction energy is considered. 
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1 Introduction 



X 

■ In a heavy-ion fusion reaction, with increasing atomic numbers of the target 

and projectile, it becomes more difficult to produce a compound nucleus owing 
to the strong Coulomb repulsion and strong dissipation forces. In the super- 
heavy mass region, this difficulty is more pronounced. Although the mecha- 
nism of the fusion-fission reaction in the heavy-mass region is not clear, gen- 
erally we recognize 

the existence of fusion hindrance ai a. 

The dynamical aspects of the fusion-fission mechanism have been investigated 
as a time evolution of the nuclear shape with friction. In the 1980s, a mean 
trajectory calculation in the nuclear deformation space using a stochastic equa- 
tion^known as the extra-push model, was performed by Swiatecki and cowork- 
ers In the superheavy- mass region, however, even if substantial extra 
push energy is supplied, it is very difficult for the mean trajectory to reach 
the spherical region owing to a s_trong dissipation. It is necessary to apply 
fluctuation-dissipation dynamics 
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The simplified calculation that included the one-dimensional diffusion model 
was applied to cold fusion reactions With one adjustable parameter, the 
results showed a good agreement with many experimental data, not only max- 
imum cross sections, but also the optimum energies and widths of the exci- 
tations functions. In the present study, we concentrate to discuss the fusion 
process using the more realistic model. 



Fusion hindrance is mainly caused by the macroscopic properties of the collid- 
ing partner (jQ). On the contrary, fusion is enhanced by the deformation and 
shell structure of the nuclei, which we call fusion enhancement (|1;E3; U; 



In experiments, especially in the case of a cold fusion reaction, these advan- 
tages are used to synthesize superheavy nuclei (flii 14). To understand the 
fusion mechanism clearly, it is better to treat separately fusion hindrance and 
fusion enhancement, that is to say, the macroscopic and microscopic aspects, 
respectively. 



In our previous study (|80, we focused, in particular, on fusion hindrance, in the 
superheavy mass region. We have presented the origin of the fusion hindrance 
systematically by conducting a trajectory calculation in three-dimensional co- 
ordinate space. To see the fundamental mechanism of the fusion hindrance, 
we employed the potential energy of the liquid drop model, and analyzed the 
mean trajectory. From the behavior of the trajectory, we can understand the 
mechanism that the fusion probability decreases exponentially as the Z num- 
ber of the fused system increases. It is concluded that the fusion hindrance 
in a system with Z greater than 102 is explained by considering the fragment 
deformation to be an important factor. 



Considering aspects of the fusion enhancement, we focus on the influence of 
the shell effect. The nuclear structure of the projectile-target combinations is 
related to the touching probability, but it also influences the dynamics from the 



toucnmg pomt to tne compound nucleus m tne superneavy mass region 
\vk s. In our approach (H;E1;@), as parameters related to the microscopic effects 
in principle we can introduce the shell correction energy on the potential 
energy surface and the transport coefficients calculated using the microscopic 
model (1^ 12; 3)- In the former case, we can see 'cold fusion valleys' in the 
potential energy surface, which have been suggested in references (3 13) • It 
is said that these valleys lead to the enhancement of the fusion probability. 
Here, we investigate precisely how the trajectory is affected by shell correction 
energy, and how the fusion probability is changed. 



In fact, shell correction energy depends on the nuclear temperature. In a dy- 
namical process, the nuclear temperature changes owing to the dissipation 
of the kinetic energy of the relative motion. Therefore, we must discuss the 
fusion process taking into account the time evolution of the potential energy 
surface, which is controlled by nuclear friction, since the dissipation depends 
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on the nuclear friction. We clarify the effect of the temperature-dependent 
shell correction energy in the fusion process. 

For clarity, in this paper, we consider the 'dynamical process' as the whole 
process from the contact point to the spherical nucleus or re-separated one. 
It includes the fusion, quasi-fission and deep quasi-fission processes, which are 
defined in reference ((6(). 

In section 2, we briefly explain our framework for the study and the model. 
We precisely investigate the effects of cold fusion valleys in section 3. In this 
section, using the potential energy surface of the full shell correction energy, 
we discuss the mechanism of the enhancement of the fusion probability by cold 
fusion valleys. In section 4, taking into account the temperature dependence of 
shell correction energy, we discuss the fusion process. In section 5, we present 
a summary and further discussion to clarify the reaction mechanism in the 
superheavy mass region. 



2 Model 



Using the same procedure as described in reference (p) to investigate the 
dynamical process, we use the fluctuation-dissipation model and employ the 
Langevin equation. We adopt the three-dimensional nuclear deformation space 
given by two-center parameterization fi21i : i22i) . The three collective parameters 
involved in the Langevin equation are as follows: zq (distance between two 
potential centers), 6 (deformation of fragments) and a (mass asymmetry of 
the colliding nuclei); a = {Ai — A2)/{Ai + A2), where Ai and A2 denote the 
mass numbers of the target and the projectile, respectively. 

The parameter 6 is defined as 5 = 3(a — 6) /(2a + b), where a and b denote 
the half length of the axes of ellipse in the z and p directions, respectively as 
expressed in Fig. 1 in reference f)2ll ). We assume that each fragment has the 
same deformations as the first approximation. 6 is related to the deformation 
parameter P2, which is familiar with us, as 

= . (1) 



The neck parameter e is the ratio of the smoothed potential height to the 
original one where two harmonic oscillat or p otentials cross each other. It is 
defined in the same manner as reference (j2lh . With e < 1, the surface of two 
fragments shows the smooth curve at the connecting point of them. On the 
other hand, in the case of e = 1, the two fragments are connected with a sharp 
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point like a top of cone. In the present calculation, e is fixed to be 1.0, so 
as to retain the contact-like configuration more realistically for two-nucleus 
collision. 



The multidimensional Langevin equation is given as 



dqi f _i\ 
dpi dV 1 d 



where a summation over repeated indices is assumed, qi denotes the deforma- 
tion coordinate specified by zq, 6 and a. Pi is the conjugate momentum of g,. 
V is the potential energy, and rriij and 'jij are the shape- dependent collective 
inertia parameter and dissipation tensor, respectively. A hydrodynamical in- 
ertia tensor is adopted in the Werner- Wheeler approximation for the velocity 
field, and the wall-and-window one-body dissipation is adopted for the dissi- 
pation tensor (|23[ 0; liih . The normalized random force Ri{t) is assumed to 
be white noise, i.e., {Ri{t))=0 and {Ri(ti)Rj(t2)) = 26ij6(ti—t2). The strength 
of random force gij is given by •jijT = J2k9ij9jk, where T is the temperature 
of the compound nucleus calculated from the intrinsic energy of the composite 
system. The potential energy is defined as 

V{q, /, T) = VUq) + ^^^77^ + Vsniq, T), (3) 

2/(g) 

VLD{q) = Es{q)+Ec{q), (4) 

v^5H(g,r) = E,%(?)<f(r), (5) 



where I{q) is the moment of inertia of a rigid body at deformation q. Vlb and 
VsH are the potential energy of the finite-range liquid drop model and the shell 
correction energy taking into account the temperature dependence. E^j^^^ is 
the shell correction energy at T = 0. The temperature dependent factor $(T) 
is discussed in the section 4.1. Es and Ec denote a generalized surface energy 
(|26l ) and Coulomb energy, respectively. The centrifugal energy arising from 
the angular momentum I of the rigid body is also considered. The detail is 
explained in reference (0). The intrinsic energy of the composite system Eint 
is calculated for each trajectory as 

E^nt = E*-\ {m-'),^p,p, - ViqJ,T), (6) 



where E* denotes the excitation energy of the compound nucleus. E* is given 
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by E* = Ecm — Q, where Q and Ecm denote the Q— value of the reaction and 
the incident energy in the center-of-mass frame, respectively. At t = 0, each 
trajectory starts from the contact point which is defined as Rtouch = Rp + Rt, 
where Rp and Rt are the radii of projectile and target, respectively. 



As discussed in reference (j60, we assumed that the kinetic energy of the rel- 
ative motion does not dissipate during the approaching process. However, to 
investigate the effect of the energy dissipation during the approaching pro- 
cess, we introduce a parameter EfUss which is the intrinsic energy converted 
from the relative kinetic energy before both nuclei touch. The discussion on 
this issue is given in section 4.3. At the contact point, the initial velocity is 
directed only in the —z direction, and the components of the initial velocities 



along the 6 and a directions are both assumed to be zero (p; l27l ) 



3 Effects of shell correction energy in fusion process 



In our previous study (jSl), we focused on the fusion hindrance from the macro- 
scopic point of view. We employed the potential energy of the liquid drop 
model and analyzed the mean trajectory excluding the final term on the right- 
hand side of Eq. (2) and used $(T) = in Eq. (5). We presented the origin 
of the fusion hindrance systematically. 

Here, we investigate the fusion enhancement from the microscopic point of 
view. The nuclear structure of the projectile-target combinations is related to 
the touching probability, but it also influences the dynamics from the contact 
point to the compound nucleus ([llLE3). We focus on the dynamics taking into 
account the shell effect after both nuclei touch. Within our model (j^ [3; 0), 
we introduce the potential energy surface with the shell correction energy. 

Figure 1 shows the potential energy surface of the liquid drop model (a) and 
with shell correction energy (6) for '^^'^llA in the z — a space (6 = 0), which is 
calculated using the two-center shell model code lioh . When we consider 
the shell correction energy on the potential energy surface, we can see pro- 
nounced valleys, that lead to the compound nucleus. The valleys are called 
'cold fusion valleys' It is said that these valleys enhance the fusion 

probability. 



We can see such shell structure in another calculation (jSCi ). Single particle 
level diagrams in reaction ™Zn+^°^Pb were calculated by a macroscopic- 
microscopic model. It shows that the entrance-channel fragment-shell effects 
remain far inside the touching point. 

In this section, we discuss the effect of shell correction energy on the dynamical 
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Fig. 1. Potential energy of the liquid drop model (a) and with shell correction energy 
(6) for ^^^114 in the z — a space (6 = 0). The calculation is done using the two-center 
sheh model code HUl). 

process using trajectory calculation. 



3. 1 Effect of cold fusion valleys on fusion process 



As our first attempt, we employ the potential energy of the liquid drop with the 
full shell correction energy, which corresponds to the potential energy surface 
at T = MeV. This potential energy is represented by Vld + Vsh{<1,T = 0), 
or Vld + E^heii by Eq. (5). We mainly discuss head-on collision, that is to 
say, the process for / = 0, because the general mechanism of the trajectory's 
behavior does not change essentially in any of the angular momentum cases, 
although the potential landscape changes with angular momentum. 



6 




Fig. 2. Potential energy of Vld and Vld + ^sheii 114 at z = 1.6, which are 

denoted by the red and blue lines, respectively. The calculation is done using the 
two-center shell model code. The combination of the Pb-target and Pu-target are 
indicated. 

Figure 2 shows the potential energies Vld and VLo+E^heii ^^^114 aX z = 1.6, 
which the potential energies correspond to a location near the contact point. 
The red and blue lines denote Vld and Vld + E^heiiy respectively. The shell 
correction energy in the reaction system where the target corresponds to Pb 
is significant, and is indicated in Fig. 2. Actually, a Pb target is chosen in a 
cold fusion reaction, in order to suppress the excitation energy of a compound 
nucleus fiSli) . 

To investigate the effects of the cold fusion valleys, we compare the fusion 
probabilities at different injection points, where the shell correction energy 
has very large negative and positive values, such as a = 0.46 and 0.24 in 
Fig. 2, respectively. 

The former corresponds to the reaction ^^Ge + ^°^Pb, which is a cold fusion 
reaction, and the latter corresponds to ^°^Ru + ^'^^Yb. Figure 3(a) shows 
sample trajectories that are projected onto the z — a {6 = 0) plane of the 
potential energy surface for ^^"^114. The contact point is marked by (+). The 
light yellow line denotes the trajectory with the starting point ao = 0.46 at 
the incident energy corresponding to the excitation energy of the compound 
nucleus ii^*=20 MeV. This trajectory appears to move along the valley till 
z ~ 1.0. Then it approaches the z ~ region. 

On the other hand, the trajectory with the starting point ao = 0.24 moves till 
2 ~ 1.0 without changing the mass asymmetry parameter a, which is denoted 
by the white line. It seems that the trajectory overcomes the mountain located 
aX z = 1.4 and a = 0.24. This trajectory does not move along the bottom of 
the valley. We can say that this trajectory is not influenced by the valley till 
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Fig. 3. Sample trajectories projected onto the z — a plane at 5 = 0.0 (a) and the 
z-b plane at a = 0.0 (b) of Vld + E^y^di ^^*^114. The light yellow and white 
lines denote the trajectories which start at a = 0.46 and 0.24 at E* = 20 MeV. 
Symbols are given in the text. 

z ~ 1.0. Then it turns in the +z direction, which corresponds to the fission 
region. 

We project the trajectories of Fig. 3(a) onto the z — 5 plane at a = 0, as shown 
in Fig. 3(b). The trajectory with = 0.46 can approach the compact shape 
region, but it is not trapped by the pocket at the ground state. The trajectory 
with ao = 0.24 moves rapidly in the +5 direction and enters the fission region. 
We can say that the former is a deep quasi-fission process (DQF) and the 
latter is a quasi-fission process (QF) (p;lsl)- 

To understand why these two trajectories take different paths, we precisely 
analyze the trajectory's behavior with the potential energy surface near the 
contact point. The potential landscape near the contact point is very impor- 
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Fig. 4. The trajectories in Fig. 3(b) are projected on to the z — b plane near the 
contact point, (a) a = 0.46 and (b) a = 0.24. The dashed hne denotes the ridge hne 
of the potential energy surface. 

tant for deciding the destiny of the trajectory ([isi). Figure 4(a) shows the 
trajectory with = 0.46, which is projected onto the z — S plane (a = 0.46). 
The trajectory is the same as that in Fig. 3. The dashed line denotes the ridge 
line. The contact point is far from the ridge line, but it does not mean the 
trajectory moves into the fission region immediately. On the potential energy 
surface near the contact point, although we can see the steep slope in the +S 
direction, there are no significant barriers in the —z direction from the contact 
point. Therefore, first the trajectory moves in the —z direction while main- 
taining 5 ~ 0. Owning to the shell correction energy, a flat potential region 
exists around z ~ 0.8 and 6 ~ 0.25. This fiat region prevents the trajectory 
from moving in the +S direction rapidly. Finally, the trajectory can enter the 
region near z ~ 0. 

Figure 4(b) shows the trajectory with ao = 0.24, which is projected onto the 
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f* (MeV) 



Fig. 5. Fusion probabilities for the initial values oq = 0.46 and oq = 0.24, which are 
denoted by the red and blue lines, respectively. 

z — 6 plane {a = 0.24). The contact point is located near the ridge line, but a 
substantial barrier exists in the —z direction from the contact point. To avoid 
the barrier, the trajectory moves in the +S direction rapidly. 

Figure 5 shows the fusion probabilities for = 0.46 and = 0.24, which are 
denoted by the red and blue lines, respectively. The potential energy at the 
contact point for a = 0.46 is lower than that for a = 0.24, as shown in Fig. 2. 
Therefore, at low excitation energy, the former fusion probability is larger than 
the latter one because in the former case the available kinetic energy at the 
contact point is larger than that in the latter case. The arrow indicates the 
Coulomb barrier (|32l). At high excitation energy, the fusion probabilities of 
these two cases are not markedly different because the kinetic energy at the 
contact point is rather large in both cases. 

3.2 Evolution of cold fusion valleys on dynamical process 

Usually, cold fusion valleys are discussed using the potential energy surface on 
the z — a plane with 5 = 0, as in Fig. 1(b). However, in the dynamical process, 
the trajectory moves in the large +6 direction. Therefore, we should discuss 
the cold fusion valleys with the dynamical evolution of the 6 parameter. For 
the 6 = case, we can see that the cold fusion valleys are remarkable. With 
changing 6 in the dynamical process, the cold fusion valleys also change. While 
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Fig. 6. Potential energy surface Vld + ^sheii z — a space at (a) 6 = 0.2 and 

(b) 6 = 0.4 for ^^^114 with a trajectory. 

the trajectory moves from z ~ 1.25 to 0.5, 6 changes to 0.2. Figure 6(a) shows 
the potential energy surface on the z — a plane at 5 = 0.2 with the trajectory. 
The position and depth of the cold fusion valleys slightly change. After moving 
along the valley, the trajectory enters the region z < 0.5. In this region, S is 
more than 0.4. Figure 6(b) shows the potential energy at 5 = 0.4 with the 
trajectory. We can see that the cold fusion valleys have already disappeared. 

3. 3 Roles of shell effects on fusion process 

Owing to the shell correction energy, we discuss another important effect to 
enhance the fusion probability. In the discussion on fusion hindrance (j^, we 
indicated that the turning point is very important. Figure 7(a) shows the 
potential energy surface of Vld for Z=102 at the turning point on the z — 6 
plane, that is to say, a corresponds to the value at the turning point. The mean 
trajectory is denoted by the white line. The red line denotes the trajectory 
taking into account the fluctuation. At the turning point {z ~ 0.3, ^ ~ 0.5), 
the trajectory moves in the fission direction owing to the potential landscape 

a. 

However, when we take into account the shell correction energy shown in 
Fig. 7(c), a temporary pocket appears at the turning point (indicated by A), 
which we call the second pocket (Il5l ). Figure 7(d) shows the potential energy 
at a = 0. The pocket (indicated by B) corresponds to the ground state. The 
trajectory is trapped in the second pocket A at the turning point, and is 
blocked from going to the fission area ([isl). While the trajectory remains in 
pocket A, the mass asymmetry is relaxed, and the large pocket B appears at 
a ~ 0. Then, the trajectory moves to pocket B. Pocket A helps the trajectory 
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Fig. 7. Sample trajectories projected onto the z — 5 plane at a which corresponds 
to the turning point, (a) Vld case and (c) Vld + ^sheii ^^^^ ^^^No. Vld + E^heii 
with a = 0.0 is shown in (d). The evolution of nuclear shape is presented in (b). 



enter the region corresponding to the compound nucleus. 



Figure 7(b) shows the time evolution of nuclear shape. The horizontal axis 
denotes 6. After the contact point, the trajectory moves in the large 6 direction 
rapidly, then it is trapped in the second pocket A. At this moment, the neck of 
nuclear shape is rather large owing to the large S, so that the mass asymmetry 
is relaxed easily while the trajectory remains in pocket A. As a approaches 
zero, the trajectory moves into pocket B. It corresponds to the compound 
nucleus. 
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Fig. 8. Sample trajectory projected onto the z — 6 plane at a = 0.42 of Vld + l^sheii 
for (a) ^''^114 (deformed nucleus) and (b) ^^^114 (spherical one). The trajectory 
starts at a = 0.67 and E* = 40 MeV. Symbols are given in the text. 

3.4 Difference between spherical and deformed compound nuclei 

In the previous section, we pointed out the role of the second pocket in the 
fusion process. However, the second pocket does not always enhance the fu- 
sion probability. As regards the position of the second pocket, we discuss the 
difficulty of synthesizing superheavy elements. 

In general, the position at the ground state of the compound nucleus is decided 
by the shell correction energy. As discussed above, in the reaction ^^Ca-f ^°^Pb, 
the ground state of the compound nucleus ^^^No is located at 5 ~ 0.2, which 
is deformed. Therefore, the second pocket A is located at a 5 with a more 
deformed shape {S ~ 0.4), as shown in Fig. 7(c), which simply corresponds to 
the turning point. In this case, pocket A assists the trajectory in entering the 
region corresponding to the compound nucleus. 

However, in the superheavy mass region {Z > 110), the nuclear shape at 
the ground state is spherical, which has a large shell correction energy. For 
example, in the reaction ^^Csi+ ^^^Pu (a = 0.67), the ground state of the 
compound nucleus 2^^114 is located at 5 = 0. In this case, the second pocket is 
located at (5 ~ 0.2, which does not correspond to the turning point {S ~ 0.4). 
Therefore, many trajectories return in the +z direction. This is one of the 
difficulties in making spherical compound nuclei such as ^^^114. 
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To show the role of the second pocket, we compare the trajectories forming the 
compound nuclei ^^^114 (deformed) and ^^^114 (spherical). Figures 8(a) and 
(b) show sample trajectories, which are projected onto the z — 5 plane at a = 
0.42 of Vld + E^f^^ii, for ^™114 and ^^^114, respectively. Both trajectories start 
at a = 0.67 with the incident energy corresponding to the excitation energy 
of the compound nucleus, ii^*=40 MeV. In this case, a ~ 0.42 corresponds to 
the turning point. 

In Fig. 8(a), a small pocket exists at z ~ 0.25 and 6 ~ 0.5, which is the 
second pocket (indicated by A). The trajectory is trapped by pocket A and 
is prevented from moving into the reseparation region. Then, it moves to the 
main pocket, which is located at around z ~ 0, 5 ~ 0.2 and a ~ 0. The 
shell correction energy at the ground state is about 5.0 MeV; therefore, the 
trajectory easily escapes from the main pocket. In Fig. 8(b), which corresponds 
to the reaction ^^Ca+^'^'^Pu, the second pocket does not exist at around z ~ 
0.25 and 6 ~ 0.5 because the compound nucleus ^^^114 is spherical. We can 
see that it is difficult for the trajectory to enter the spherical region. The 
trajectory returns to the fission region. 

Figure 9 shows the fusion probabilities for forming the compound nuclei ^^^114 
and ^''^114, which are denoted by the blue and red fines, respectively. We use 
the initial value of a = 0.67 for both cases. The fusion probability of ^''^114 
is larger than that of ^^^114, owing to the second pocket, which is located at 
the deformed shape. 



4 Effect of temperature-dependent shell correction energy 



In the previous section, we used the full shell correction energy for the potential 
energy surface. However, the shell correction energy actually depends on the 
nuclear temperature T. For the dynamical calculation for the fusion process, 
the temperature dependence of shell correction energy has not been considered 
sufficiently in many publications. 

In this section, we discuss the temperature dependence of the shell correction 
energy and how the fusion process is affected. Considerable effort has been 
made to investigate the temperature dependence of the level density parameter 



33l : |3J) and applied to the calculation of the statistical model for the fission 



process fl35f l. The temperature dependence of the potential energy surface has 



been investigated in references (|36l : |37|). Here, we discuss the effect of this 



dependence on the fusion process on the basis of dynamical calculations. 
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Fig. 9. Fusion probabilities for forming the compound nuclei ^^^114 (spherical) and 
^^^114 (deformed), which are denoted by the blue and red lines, respectively. We 
use the initial value of q = 0.67 for both cases. 

J^.l Temperature dependence of shell correction energy 



The temperature dependence of shell correction energy can be extracted from 
the free energy calculated using single particle energies, which is obtained from 
the two-center shell model code fj28i : l29i ). The procedure is given precisely in 
reference (|38h. The temperature dependence of shell correction energy calcu- 
lated by this method is represented by EsheuilyT). 

The ratio of shell correction energies, EgheiiiQ, T)/ Egheiiiq, T = 0), in the spher- 
ical region for ^^^114 is shown in Fig. 10(a) by the solid line. As discussed in 
reference (jsih. our result is consistent with the parameterization of the factor 
$(T), 



$(r) = exp 



aT2 



(7) 



following the work by Ignatyuk et al. ("s^, where a denotes the level density 
parameter of the constant value a = A/8. The shell damping energy is cho- 
sen as 20 MeV, which is given by Ignatyuk et al. (fj^. This parameterization 
is denoted by the open circles in Fig. 10(a), which shows a tendency similar 
to that of the solid line. 
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Fig. 10. Calculated temperature dependence of shell correction energy, (a) For '^^'^114: 
at g = (the spherical shape), compared with Ignatyuk's parameterization with 
a = A/8 and Ed = 20 MeV. (b) For ^^Ge+^^^Pb reaction at the contact point. The 
symbols and lines are explained in the text. 

In the fusion process, it is important to determine EsheuiliT) near the con- 
tact point of the colhding partner. The potential landscape near the contact 
point mainly decides the destiny of the trajectory in the dynamical process. 
Using the same procedure, we calculate the temperature dependence of shell 
correction energy at the contact point Esh^uiq = cont,T) for the reaction 
76Qg^208p|3_ In Fig. 10(b), the blue line denotes Esheii{q = cont,T)/ Esheiiiq = 
cont,T = 0). This tendency is different from that of the spherical shape case, 
which is denoted by the dashed line. In Eq. (7), when we take into account the 
nuclear deformation dependence of the level density parameter ax.sXl) (39), 
our result can be fitted using E^ = 20 MeV, denoted by the open circles. The 
result shows the similar tendency of the blue line. 



4-2 Trajectory behavior on temperature- dependent potential energy surface 



In the fission process, the time evolution of the potential energy surface affects 
the dynamical process considerably. Owing to the neutron emission from the 
compound nucleus, the temperature of the nucleus decreases, and the shell 
correction energy recovers. The time scale of the variation of the potential 
energy surface affects the dynamics of the fission process (jsih . 
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Fig. 11. Time evolution of <I>(T) for a trajectory in the reaction ''^Ge+^'^^Pb at 
E* = 40MeV with E^iss = MeV. 

Here, we focus on the effect of the temperature-dependent shell correction 
energy on the fusion process. In the dynamical calculation, we analyze the 
trajectory on the temperature-dependent potential energy surface, which is 
concerned with the intrinsic energy Eint of the system. 

Within our model, the time evolution of E^nt is controlled by the friction 
tensor 7. Also, E^t is connected with Ediss at the contact point (ji?! ). Here, 
we clarify the relationship among these variables and the mechanism of the 
fusion process. 

We start the dynamical calculation from the contact point using the Langevin 
equation. Ediss is treated as a parameter of the initial condition. Discussions 
on this issue have been given in references (4^ 13). In fact, E^iss is related to 



the dynamics in the approaching process. Considering EfUss as a parameter, 
we investigate the effect of E^iss on the trajectory in the process after the both 
nuclei touch. Then, in the next section, we discuss the effect of the friction 
tensor 7. 

Before we discuss the effect of the different Ediss's on the dynamical process, we 
analyze the trajectory under the influence of the temperature-dependent shell 
correction energy, using the same initial condition as Ediss = MeV in the 
previous calculation. We employ the potential energy surface Vl£)+Vsh{T = 0) 
and Vld + Vsh{T) in the reaction ^Ge+^o^Pb. 

When we use Vld + VsHiT) as the potential energy surface, as the intrinsic 
energy that is transformed from kinetic energy increases, some part of the shell 
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Fig. 12. Sample trajectory on Vld + Vsh{T = 0) and Vld + Vsh(T) in the reaction 
^^Ge+^''^Pb at E* = 40 MeV, which are denoted by blue and red lines, respectively. 
The trajectories are projected onto (a) the z — a plane and (b) the z — 5 plane, 
respectively. Symbols are given in the text. 

correction energy disappears. Figure 11 shows the time evolution of $(T) for 
a trajectory at E* = 40 MeV. At the contact point, where t = in the 
calculation, owing to the assumption Ediss = MeV, which corresponds to 
Eint = MeV, the full shell correction energy is indicated. Then, till t ~ 
1.0 X 10^^^ sec, more than 50% of the shell correction energy disappears owing 
to the strong friction tensor. During this extremely short reaction time, the 
kinetic en ergy dissipates into the intrinsic energy, which has been discussed in 
reference (|l5ri. 

As consequence, such variation of the potential energy surface gives the in- 
fluence on the trajectory's behavior. Figure 12 shows sample trajectories on 
Vld + Vsh{T = 0) and Vld + Vsh(T), which are denoted by the blue and 
red lines. The trajectories are projected onto the z — a plane and z — 6 plane, 
which are shown in (a) and (b), respectively. The (+) denotes the point of 
contact in the system. Both trajectories in Fig. 12 are calculated using the 
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Fig. 13. Fusion probabilities Pcn in the reaction ^^Ge+^'^^Pb for each potential en- 
ergy with Ediss = MeV. The blue and red lines denote Pcn with Vld+Vsh{T = 0) 
and Vld + VsHiT), respectively. Pcn with Vld is represented by the black line. 



same random numbers. After t = 7.9 x 10~^^ sec, which is indicated by the 
black arrow, trajectories separate along different paths owing to the different 
potential energy surfaces. The trajectory with Vld + Vsh(T = 0) takes the 
fusion process. On the other hand, the trajectory with Vld + VsniT) takes 
the deep quasi-fission process. This is brought about by the disappearance of 
the second pocket at around z ~ 0.2, 6 ~ 0.4 and a ~ 0.3. 

Figure 13 shows the fusion probabilities Pcn in the reaction ^^Ge+^°^Pb for 
each potential energy with E^iss = MeV. The blue and red lines denote 
Pcn with Vld + Vsh(T = 0) and Vld + Vsh(T), respectively. Pcn with Vld 
is represented by the black line. With increasing incident energy, Pcn with 
Vld + Vsh(T) approaches that with Vld owing to the disappearance of shell 
correction energy at high temperature. 

The mass distributions of fission fragments in the reaction ^Ge+^^'^Pb at 
E* = 40 MeV are shown in Fig. 14. The blue and red lines denote the cases 
of Vld + Vsh{T = 0) and Vld + VsniT) with E^iss = MeV, respectively. 
In the latter case, some part of the shell correction energy disappears and the 
number of mass symmetric fission events increases in comparison with that in 
the Vld + Vsh{T = 0) case. 
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Fig. 14. Mass distribution of fission fragments in the reaction Ge-|- Pb at 
E* = 40 MeV. The blue and red hues denote the cases of Vld + Vsh(T = 0) 
and Vld + Vsh{T) with Ediss = 0%, respectively. The gray shadow denotes the 
case of Vld + Vsh{T) with E^iss = 100%. 

4-3 Ediss dependence on dynamical process 



Next, we focus on the dependence on E^iss- As discussed in references f)27l: l40[). 
the fusion probability and mass distribution of fission fragments are affected by 
Ediss- When we use the condition Ediss = MeV (or Ediss = 0%), we assume 
that the velocity of the trajectory is in the —z direction at the contact point, 
so that there exists trajectories that enter the fusion region. However, when 
we take the condition in which all the kinetic energy has already dissipated 
at the contact point, which is represented by Ediss = 100%, all trajectories 
move down along the steep slope in the fission direction (the +z direction) 
immediately. As a result, the fusion probability is equal to zero. We discuss 
this point in connection with the potential energy surface near the contact 
point in section 4.5. 

First, we investigate the evolution of $(T) related to Ediss- When the kinetic 
energy dissipates during the approaching process, some part of the shell cor- 
rection energy has already disappeared at the contact point. Figure 15 shows 
the time evolution of $(T) for a trajectory at E* = 40 MeV with various 
Ediss- At t = sec, $(T) strongly depends on Ediss- In this case, the potential 
energy surface at the contact point for each Ediss is shown in Fig. 16. Owing 
to the strong friction, 60% of the shell correction energy is rapidly lost by 
7.5 X 10~^^ sec. After this time, the trajectory of each Ediss nioves on almost 
the same potential energy surface. 
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Fig. 15. Time evolution of $(T) for a trajectory in the reaction ^^Ge+^ospb at 
E* = 40 MeV for various Ediss cases. 



In Fig. 15, the gray hne denotes $(T) for Ediss = 100%, which corresponds 
to Ediss = 34 MeV in this case. In fact, in the surface friction model, it shows 



that all the kinetic energy dissipates during the approaching process ()41l: |4 



Under this condition, the gray line shows that 50% of the shell correction 
energy has already disappeared at the contact point. The trajectory moves 
down the steep slope rapidly and arrives at the fission region at 8.0 x 10~^^ 
sec. The mass distribution of the fission fragments in this reaction is denoted 
by the gray shadow in Fig. 14. As discussed in reference (|4ol). it shows the two 
sharp peaks corresponding to target-like and projectile-like fragments. 



The fusion probabilities for the various Ediss cases are shown in Fig. 17. We 
can roughly say that the fusion probability of each Ediss shifts by as much 
as Ediss MeV from the black line in the Ediss = MeV case. For example, 
in Fig. 17, we connect the points that have the same value of E* — Ediss at 
the contact point using the red dashed line. On this dashed line, the fusion 
probability with the large Ediss increases owing to the thermal fluctuation of 
the collective motion. However, we can see the systematic tendency and the 
relationship between the fusion probability and Ediss- Even though we treat 
Ediss as a parameter in our model, the mechanism of the dynamical process 
after the contact point is clarified. 
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Fig. 16. Potential energy surface Vld + VsniT) at the contact point in the reaction 
^'^Ge+^ospb at E* = 40 MeV for various E^iss cases. 
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Fig. 17. Fusion probabihties Pcn in the reaction '^'^Ge+^O'^Pb with Vld + VsniT) 
for various Ediss cases. 
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Fig. 18. Time evolution of $(T) for the different friction tensors 7 in the reaction 
76Qg_^208p|3 at E* = 40 MeV with Ediss = MeV. 

4-4 Friction tensor dependence on dynamical process 



In the dynamical process, the dissipation speed is controlled by the friction 
tensor 7. That is to say, the time evolution of Eint depends on 7. The ambiguity 
of the friction tensor derived from the macroscopic model has been discussed 



very often ()43t l44l : 1451 ) . Moreover, as for the effect of nuclear structure, we 



should consider the microscopic transport coefficients, that are calculated by 



the linear response theory ()16l : llTt Il8^ . The friction tensor that calculated 



by the microscopic model depends on the temperature and is smaller than 
that calculated by the macroscopic model at low temperature. Here, using the 
modified strength of the one-body friction tensor, we investigate the effect of 
the friction tensor on the dynamical process, in terms of the variational speed 
of the potential energy surface. 

Figure 18 shows the time evolution of $(T) for the different friction tensors 
in the reaction ™Ge+^°^Pb at E* = 40 MeV with E^iss = MeV. The factors 
multiplied with the friction tensor are 0.1, 0.5, 1.0, 1.5 and 2.0, and the cor- 
responding results are shown. Except in the case of the multiplied factor 0.1, 
more than 40% of the shell correction energy disappears till 3.0 x 10^^^ sec 
for all cases. This time is rather short in comparison with the reaction time 
in the whole dynamical process. After t = 5.0 x 10~^^ sec, in the smaller fric- 
tion case, $(T) decreases faster than in the larger friction case. In the smaller 
friction case, after the trajectory overcomes the fusion barrier rapidly, the tra- 
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Fig. 19. Fusion probabilities Pcn for the different friction tensor 7 in the reaction 
76Qg^208p^j with Ediss = MeV. Vld + Vsh{T) is employed. 



jectory enters the region where the potential energy is lower than that near the 
contact point. It is possible that the variation in the potential energy is trans- 
formed into the intrinsic energy by Eq. (6). In this way, after the trajectory 
overcomes the fusion barrier at high speed, the intrinsic energy increases and 
the shell correction energy disappears rapidly. The intrinsic energy increases, 
not only owing to the dissipation from kinetic energy, but also owing to the 
transformation of the variation from the potential energy. 



Figure 19 shows the fusion probabilities Pcn for each friction tensor for the 
reaction ''^Ge+^^^Pb with Ediss = MeV. Pqn increases with decreasing mag- 
nitude of the friction tensor. The friction tensor affects the dynamical process 
considerably (|27l). 



We investigate the mass distribution of fission fragments for each friction ten- 
sor for the reaction ^^Ge+^ospb at E* = 40 MeV with E^iss = MeV, which is 
shown in Fig. 20. In the small friction case, mass symmetric fission events are 
dominant, because the trajectory can reach the small z region. On the other 
hand, mass asymmetric fission events are dominant in the large friction case. 
The smallest value of z that the trajectory can reach is related to the mass 
distribution of fission fragments ([ish. 
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Fig. 20. Mass distribution of fission fragments for each friction tensor 7 in the 
reaction '^^Ge+^ospb at E* = 40 MeV with Ediss = MeV 

4-5 Effect of modified potential energy surface 



As we mentioned in section 4.3, when all the kinetic energy dissipates during 
the approaching process {Ediss = 100%), all trajectories move down the steep 
slope in the fission direction and the mass distribution of fission fragments 
has two sharp peaks, as shown in Fig. 14. In connection with this problem, we 
should precisely investigate the landscape of the potential energy surface near 
the contact point. The modification of the potential energy near the contact 
point is performed continuously (jiih . 

As a problem in our two-center shell model, since we employ the mass asymme- 
try parameter a, the charge distribution of the nuclei follows a. In the fusion 
process, this is not the proper assumption before the contact point, which we 
call the two-body region. At least, in the two-body region, we should calcu- 
late the potential energy using realistic distributions for the proton charge 
az and mass um of the projectile and target. By performing the analysis in 
this manner in the two-body region, we obtain the modified potential energy 
surface by employing the proper smoothing between the two-body region and 
the one-body region at the contact point, which is denoted by the blue line in 
Fig. 21. It shows the ^^Ca+^^'^Pu reaction case, where both colliding nuclei are 
assumed to have a spherical shape. The red line denotes the potential energy 
using the same distributions of a for the charge and mass of both nuclei, even 
in the two-body region, that were used in the previous calculation. The open 
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Fig. 21. Potential energy in the reaction ^Ca+ Pu calculated using the two-center 
shell model. The red line denotes the potential energy using om = c^z for all region. 
The blue line denotes the modified potential V"^°'^ using aM and az of the projectile 
and target in the two-body region. The potential energy calculated using the Bass 
model is denoted by the green line. Symbols denote in the text. 

circles denote the contact point of the colliding partner. The Coulomb barrier 
of the modified potential energy agrees with that calculated using the Bass 
model, which is denoted by the green line (H^)- 

Using the modified potential energy surface, which is represented by V'^°'^, 
we calculate the fusion probability and mass distribution of fission fragments 
using the same procedure as in the previous section. We start the Langevin 
calculation at the contact point. We also take into account the temperature 
dependence of shell correction energy. 

In the case of Ediss = MeV, the fusion probability with the modified potential 
ymod ^Qgg j-^Q^ change considerably compared with the previous results. Since 
the velocity of the trajectory is in the —z direction, it is easy to enter the 
small z region, and it is not affected by the potential energy surface in the 
two-body region. In the case of E^iss = MeV, Fig. 22 shows Pcn using the 
previous potential and the modified potential, which are denoted by the red 
and blue lines, respectively. 

On the other hand, when all the kinetic energy dissipates during the approach- 
ing process, even though the fusion probability with the previous potential is 
zero for any incident energy, it markedly increases in the case of the modified 
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Fig. 22. Fusion probabilities Pcn in the reaction ''^Ca+^^^Pu for each potential. In 
the case of Ediss = MeV, Pcn using the previous potential energy V is denoted by 
the red line. Pcn using the modified potential energy 

ymod £qj. £j^g different E^iss 

are denoted by each line. 



potential V"^""^, which is denoted by the green line in Fig. 22. In the V"^°'^ 
case, the Coulomb barrier located at z = 1.63 prevents the trajectory from 
moving into the fission region immediately, and the trajectory at the contact 
point moves down along the potential slope in the —z direction automatically. 
In this case, the trajectory has the momentum to go to the fusion region. In 
Fig. 22, we also show Pcn in the Ediss = 50% and 75 % cases with V"^"'^. Pcn 
decreases with increasing Ediss- 

Figure 23 shows the mass distribution of fission fragments with the modi- 
fied potential energy surface, in the reaction ^^Ca+^^^Pu at E* = 36.5 MeV. 
When all the kinetic energy dissipates during the approaching process, the 
results with the previous and modified potential energy surfaces are denoted 
by the red shadow and green lines, respectively. Whereas the result with the 
previous potential shows two sharp peaks, that with the modified potential 
is distributed around the wide mass region. The tendency of the latter shows 
a rather good agreement with the experimental data, which correspond to 



fusion-fission and quasi-fission events, as denoted by the black line (1471 ). Here, 
the experimental data are normalized to agree with the results in the V"^°'^ 
case with Ediss = 100% at the mass symmetric point. Also, we plot the results 
of different Ediss cases in Fig. 23. 

Even though this result is a preliminary one, it shows the importance of the 
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Fig. 23. Mass distributions of fission fragments with the modified potential energy 
ymod ^Yie reaction ^^Ca+^^^Pu at E* = 36.5 MeV for the different Ediss- The 
red shadow denotes the mass distributions using the previous potential energy with 
Ediss = MeV. The experimental data is denoted by black line (ji^). 

potential energy surface near the contact point and Ediss- We should inves- 
tigate more precisely the effect of the approaching process with the modified 
potential energy. 



5 Summary 

The dynamical process in the superheavy mass region was studied system- 
atically by trajectory calculation in three-dimensional coordinate space. To 
understand the fusion mechanism clearly, we treated the dynamical process as 
a fusion hindrance and as a fusion enhancement, separately. 

On the basis of our previous results we discussed the fusion enhancement. 
We introduced shell correction energy on the potential energy surface, and 
investigated how trajectory is affected by shell correction energy. 

The fusion probability increases along the cold fusion valleys at low incident 
energy. This is caused by the potential energy surface near the contact point 
on the z — 6 plane. With shell correction energy, the second pocket appears at 
the location with a deformed shape. When the position of the second pocket 
corresponds to the turning point, it enhances the fusion probability. 
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The temperature dependence of shell correction energy was also discussed. 
In the dynamical process, the evolution of $(T) is related to Ediss and the 
friction tensor 7. We analyzed the trajectory on the temperature-dependent 
potential energy surface, and clarified the effects of E^iss and 7. The fusion 
probability and mass distribution of fission fragments were calculated. 

This study is the first systematic attempt to investigate the role of shell cor- 
rection energy in the dynamical process. Shell correction energy plays a very 
important role in the fusion process and enhances the fusion probability, when 
the colliding partner has a strong shell structure. Although we treated E^iss 
as a parameter, after the both nuclei touch, the mechanism of the dynamical 
process and the relationship among the parameters are clarified. 

As a topic of further study, the potential energy surface should be improved 
near the contact point and we should investigate the dynamical process more 
accurately (|46h. We should take into account the dynamics in the approaching 
process and clarify how the kinetic energy dissipates in the approaching pro- 
cess. Also, we should directly apply the transport coefficients calculated using 
the microscopic model, for example, linear response theory 
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